The SO(5) Theory as a Critical Dynamical Symmetry 



ON 

on 

On 



Lian-Ao Wu (1) , Mike Guidry( 2 \ Yang Sun( 2 ), and Cheng-Li Wu (3) 
^Department of Physics, Jilin University, Changchun, Jilin, 130023, PRC 
^Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996 
^Department of Physics, Chung-Yuan Christian University, Taiwan, ROC 

(February 1, 2008) 

We use generalized coherent states to analyze the SO (5) theory of high-temperature superconduc- 
tivity and antiferromagnetism. We demonstrate that the SO(5) symmetry can be embedded in a 
larger algebra that allows it to be interpreted as a critical dynamical symmetry interpolating be- 
tween antiferromagnetic and superconducting phases. This dynamical interpretation suggests that 
SO (5) defines a phase with the character of a spin-glass for a significant range of doping. 
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Data for cuprate high-temperature superconductivity 
(SC) suggests D-wave singlet pairing and that SC in these 
systems is closely related to the antiferromagnetic (AF) 
insulator properties of the undoped compounds. Zhang 
et al proposed to unify AF and SC states by assem- 
bling their order parameters into a 5-dimensional vector 
and constructing an SO (5) group that rotates AF or- 
der into D-wave SC order. Recently Q, we introduced a 
U (4) model of high-temperature SC and AF order hav- 
ing SO (5) as a subgroup. In this paper we use coherent 
states H to analyze this subgroup and provide a dynam- 
ical interpretation of 50(5) and the order parameters 
that it rotates. The resulting energy surfaces suggest 
that SO (5) is a critical dynamical symmetry interpolat- 
ing between AF and SC phases that is extremely soft 
against AF fluctuations. Thus, we shall propose that the 
Zhang 50(5) theory corresponds to a spin-glass phase 
that can serve as a doorway between AF and SC order. 

Gilmore 0,7| and Perelomov |sj (see also earlier work 
by Klauder |9|), demonstrated that Glauber coherent 
states [jlOj for the electromagnetic field could be gen- 
eralized to coherent states for an arbitrary Lie group. 
These states can be analyzed in terms of their geome- 
try, which is in one-to-one correspondence with the coset 
space. However, it is often simpler to view them as 
Hartree-Fock-Bogoliubov (HFB) variational states con- 
strained by the dynamical symmetry. The symmetry- 
constrained HFB coherent state method is discussed in 
Refs. It may be viewed as a type of mean- 

field approximation to the underlying many-body prob- 
lem that is particularly useful in the present context be- 
cause it leads to easily visualized energy surfaces. This 
provides a natural connection to spontaneously broken 
symmetries and effective Lagrangian field theories. 

Let us begin with a group structure. We introduce 16 
bilinear fermion operators: 
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where c[ . creates a fermion of momentum k and spin 
projection i,j = 1 or 2 = f or I, Q = (n, n,ir) is an AF 
ordering vector, f2 is the lattice degeneracy, and following 



Refs. |2]]3| we define g(k) — sgn(cosfc x — cosk y ) with g(k+ 
Q) = -g(k) and \g(k)\ = l. 

Under commutation the operator set (|l|) closes a t/(4) 
algebra corresponding to the group structure 

D 50(4) x U(l) D SU(2) S x U{1) 
C7(4) D SU{4) D 50(5) D SU(2) S x 17(1) (2) 
D SU(2) P x SU(2) S D SU(2) S x 17(1) 

where we require subgroup chains to end in SU(2) S x Z7(l) 
representing spin and charge conservation. In Ref. 0] 
we discuss the representation structure of (Q) and show 
that the 50(4) subgroup is associated with AF and the 
SU(2) p subgroup with D-wave SC; in this paper, we jus- 
tify interpreting the 50(5) subgroup as a critical dynam- 
ical symmetry leading to a "spin glass" phase (SG). 

It is convenient to take as the generators of U(4) — » 
[7(1) x SU(A) the combinations 

0+ = Oil + 022 = E( C fc+QT Cfc T + C k+Ql°kl) 
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where 0+ generates the U(l) factor and is associated 
with charge density waves, 5 is the spin operator, Q is the 
staggered magnetization, the operators 7ft , 7? are those of 
Ref. §, the operators , D are associated with D-wave 
pairs, and M is the charge operator. 
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To facilitate comparison with the 50(5) symmetry the 
50(6) ~ SU(4) generators may be expressed as 
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with L a b = —Lb a and with commutation relations 

[L a b, L c d] = i(5 ac Lbd — SadLbc — SbcLad + <5bdL a c)- 



(5) 
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The coherent state method requires a faithful matrix 
representation of SU(4). Explicit multiplication verifies 
that the following mapping preserves the algebra of (0) : 
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where <7j are Pauli matrices, tr± = ^(c x ±«cr y ) and / is a 
unit matrix. The Casimir operator of SU(4) ~ 50(6) is 

C su(4) =7f t -7f + D t L> + 5-5+ Q-Q + Af(Af -4) (8) 

and the irreps may be labeled by the 50(6) quantum 
numbers, (01, 02, 03). We take as a Hilbert space 

|5) = \n x n y n z n d ) = (^r^""^)"^)"*! } ( 9 ) 

which is a collective subspace associated with 50(6) ir- 
reps of the form (oi, CT2, 0-3) = (§,0, 0). 

Utilizing the methods of Ref. Q, the coset space is 
5t/(4)/50(4)x£/(l), where the 50(4) subgroup is gen- 
erated by Q and 5 and C/(l) is generated by the charge 
operator M. The coherent state may be written as 

I V) - cxp(Wi2 + moq\ 2 - h. c.) I 0) = T I 0), (10) 



where the real parameters 7700 and 7710 are symmetry- 
constrained variational parameters. Since they weight 
the elementary excitation operators p\ 2 and q\ 2 in Eq. 
( pTj| ) , they represent collective state parameters for a sub- 
space truncated under the 5£/(4) symmetry. The most 
general coherent state corresponds to a 4-dimensional, 
complex, compact manifold parameterized by 8 real vari- 
ables. The reduction of the coherent state parameters to 
only two in Eq. ( p~C| ) follows from requiring spin and time 
reversal symmetries for the most general wavefunction. 

From the coset representative expressed in the 4- 
dimensional matrix representation (^), the transforma- 
tion T can be written as 
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where Yi and Y 2 are determined by the requirement that 
T be unitary and a and (3 are variational parameters re- 
lated to r] 00 and 7710 in Eq. ( |io| ) (see Ref. Q). Introducing 
v + = a + (3 and tj_ = a — 0, and u + = (1 — v^) 1 / 2 and 
u_ = (1 — tj 2 ) 1 / 2 , the matrices can be written as: 
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We may also introduce the gap parameters A + = it+w+ 
and A_ = 

One can regard Eq. (jlTJ) as a Bogoliubov type transfor- 
mation, which, through the operator T, transforms the 
bare vacuum state |0) to a quasiparticle vacuum state 
\ip). Likewise, the basic fermion operators {c^,c} can be 
transformed to quasifermion operators {a^,a} through 



T 



(12) 



Using the transformation (|l2|), one can express any one- 
body operator in the quasiparticle space jL5|. In this 
way, many results known in HFB theory can be applied 
directly in the present case. Generally, 



(13) 



with the Cs determined by the transformation properties 
of the operator O (see Appendix E of Ref. [^5[). 

The expectation value for an operator O in the coher- 
ent state representation is given by (O) = Tr(0( 22 )) and 
for a 2-body operator O^O, 



(6*6) = Tr(0^ 22) )Tr(0^ 22) ) + Tk(0t{ai)0( 



12) 
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Utilizing (|7|) and (|1 1[) , the expectation value of the par- 
ticle number in the coherent state is 



(N) = 2£l(a 2 + f3 2 ) = n{v 2 + + v 2 _). 



(15) 
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FromEq. @, the squares of a and (3 are constrained by 
the equation of a circle having radius sj n/2fi. Thus, for 
fixed particle number n we may evaluate matrix elements 
in terms of a single variational parameter, say (3, which 
may be related to standard order parameters by compar- 
ing matrix elements. For example, the z component of 
the staggered magnetization is given by 

(Q z ) = 2fi/?(n/(2fi) - (3 2 ) 1 ' 2 = ±ft(«+ - (16) 

Let us consider a simple 50 (5) Hamiltonian given by 
H$ ~ C so (5) where the Casimir for 5*0(5) is 

C so (5) =7? t - tt + S-S + M(M -3). (17) 

We evaluate the SO (5) ground-state energy surface by 
taking the expectation value of H$ in the coherent state 
(|lO|). The expectation value of the 50(5) Casimir oper- 
ator may be written as [see Eq. (J|)] 

<0 5 > = (C su{4) ) - (£>t£>) - ($ ■ 3) + (n - n)/2 (18) 

where <O su(4) } = § (f +4). From (@), ©, and (0), 

(DtU) = ±fi 2 (A+ + A_) 2 + |n(«* + (19) 

+ 5«[A+ + A 2 . + {u+v- +U-V+) 2 ], (20) 

and the energy surface follows from Eqns. (p7|)-(|2C|). 

Eq. ( |l8| ) is an expectation value of a Hamiltonian that 
can be parameterized more generally as 

H=-G \{l-p)rtD+p4-~&], (21) 

with Go > and < p < 1. Then p = \ corresponds 
to 50(5) symmetry, while the extreme values and 1 
correspond to SU (2) and SO (A) symmetries, respectively 
(see Ref. ||). Other values of p respect 5£/(4) symmetry 
but break the 50(5), 50(4), and SU{2) subgroups. In 
Fig. 1 we illustrate the classical ground-state energy (H) 
as a function of the order parameter (3 for different lattice 
occupation fractions with p = 0, ^ and 1. 

For p = [SU (2) limit], the minimum energy occurs at 
((3) = for all values of n, indicating SC order. For p = 1 
[SO (4) limit], the opposite situation occurs: (j3) = is an 
unstable point and an infinitessimal fluctuation will drive 
the system to the energy minima at finite {(3) (and thus 
finite (Qz)), causing a transition to AF order. From Fig. 
lb, the 50(5) dynamical symmetry is seen to have ex- 
tremely interesting behavior: For n near £1 (half filling), 
the system has an energy surface almost flat for broad 
ranges of (3, corresponding to large-amplitude fluctua- 
tions in AF order. We conclude that the 50(5) symme- 
try represents a phase having much of the character of a 
spin glass for a range of doping fractions. 

Under exact 50(5) symmetry, AF and SC states are 
degenerate at half filling and there is no classical barrier 



between them (the n/Cl = 1 curve of Fig. lb). To see 
clearly how the 50(5) interpolates between AF and SC 
states as particle number varies, let us perturb slightly 
away from the 50(5) limit of p — 1/2. In Fig. 2a, re- 
sults for p = 0.52 are shown. We denote the value of ((3) 
minimizing (H) as (/3)o- For small values of n the stable 
point is (/3)o = 0. This corresponds to SC order. If n 
is near fl (half filling), (/3)q ~ ±0.5; this corresponds to 
AF order. For intermediate values of n there is a sub- 
stantial region in which the system has an energy surface 
almost flat for broad ranges of (/?) , implying the presence 
of large-amplitude fluctuations in AF order. The corre- 
sponding variation of the AF order parameter ( Q z } with 
n for several values of p is summarized in Fig. 2b. 

Thus, as particle number varies the 50(5) symmetry, 
or slightly perturbed 50(5) symmetry, interpolates be- 
tween AF order at half filling and SC order at smaller 
filling, with the intermediate regime acting like a spin 
glass. Because for exact 50(5) symmetry there is no 
barrier at half filling between AF and SC states, one can 
fluctuate into the other at zero cost in energy. Only when 
50(5) is broken does the energy surface interpolate be- 
tween AF and SC order as doping is varied (compare the 
surfaces for p = 1/2 and p = 0.52). From Ref. and the 
present paper we conclude that the underdoped portion 
of the cuprate phase diagram is described by a Hamil- 
tonian that conserves 51/(4) but breaks (explicitly and 
weakly) 50(5) symmetry in a direction favoring AF order 
over SC order. Thus, we propose that the underdoped 
regime is associated with a conserved SU (4) but weakly 
broken 50(5) symmetry. Likewise, optimally doped su- 
perconductors and AF insulators are associated with the 
symmetries SU(4) D SU{2) and SU(4) D 50(4) respec- 
tively, or small perturbations around them. 

Dynamical symmetries that interpolate between other 
dynamical symmetries have been termed critical dynam- 
ical symmetries EJ. Thus, 50(5) is a critical dynam- 
ical symmetry. Such symmetries are well known in nu- 
clear structure (T^jlJ,^], and the 50(5) critical dynami- 
cal symmetry discussed here has many formal similarities 
with the SO (7) critical dynamical symmetry of the (nu- 
clear) Fermion Dynamical Symmetry Model Jl(|. The 
condition for realization of the 50(5) phase is that the 
strength of Q ■ Q equals that of D'D in the Hamilto- 
nian; This is similar to the 50(7) nuclear critical dy- 
namical symmetry, which is realized when there is an 
overall 50(8) symmetry and the pairing and quadrupole 
interaction terms are equal in the Hamiltonian [ pT[ . 

Finally, let us revisit a critical conclusion of Ref. 0] 
in light of the present discussion. In the 5C/(4) model, 
the 50(5) subgroup defines only one possible dynamical 
symmetry. The AF phase near half filling and the SC 
phase near optimal doping of the SU(4) model are as- 
sociated, not with the 50(5) subgroup, but with SU(2) 
and 50(4) subgroups, respectively. The 50(5) subgroup 
(identical to the Zhang 50(5) but embedded in a larger 
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group structure) defines a phase interpolating between 
the SU(2) and 50(4) phases and characterized by large 
AF and SC fluctuations. Thus, in the 577(4) model the 
relationship of the AF phase near half filling and the 
SC phase near optimal doping is not constrained by the 
50(5) subgroup, and SO (5) should be viewed as the ap- 
propriate symmetry for the undcrdoped compounds. 

To conclude, we have proposed separately B that high 
T c behavior of the cuprates results from a ?7(4) symmetry 
realized dynamically in terms of 3 subgroup chains. One 
of those subgroups is the 50(5) subgroup discussed ex- 
tensively by S. C. Zhang |Q. In this paper, we have used 
SU(A) coherent states to analyze the 50(5) energy sur- 
face. We interpret the Zhang 50(5) as a critical dynam- 
ical symmetry that interpolates between AF and D-wave 
SC order as doping is varied, and have noted similari- 
ties with analogous critical dynamical symmetries from 
nuclear structure physics. This permits the 50(5) sym- 
metry to be understood dynamically as a critical phase 
that for a range of doping has an energy surface extremely 
soft against AF fluctuations and therefore having much 
of the character of a spin glass. Thus, we propose that 
slightly broken 50(5) is the symmetry of the underdoped 
regime, but that AF and optimally doped SC compounds 
are described by different subgroups of SU (4) . 

L.-A. Wu was supported in part by the National Nat- 
ural Science Foundation of China. Research at Chung 
Yuan Christian University is supported by the National 
Science Council of the ROC through contracts NO. NSC 
88-2112-M033-006. 



Ginocchio, Nucl. Phys. A505, 7 (1989) 
[15] P. Ring and P. Schuck, The Nuclear Many-Body Problem, 

(Springer- Verlag, New York, 1980) 
[16] C.-L. Wu, D.H. Feng and M. W. Guidry, Adv. in Nucl. 

Phys. 21, 227 (1994) 

FIG. 1. Coherent state energy surfaces. The energy unit 
is Gof2 2 /4 (see Eq. (pH)). Numbers on curves are twice the 
lattice occupation fractions, with n/Q = 1 corresponding to 
half filling and < n/Q < 1 to finite hole doping. SO(5) 
symmetry corresponds to p — 1/2. 

FIG. 2. (a) As for Fig. 1, but for slightly perturbed SO(5). 
The dotted line indicates the classical ground state as n varies, 
(b) Variation of the AF order parameter with occupation 
number for different values of p. 
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